Ferroelectricity in Pb(Zr .5Tio.5)03 thin films: critical thickness and 180° stripe 

domains 
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The ferroelectric properties of disorder Pb(Zro.sTio.5)03 thin films are investigated with Monte 
Carlo simulations on the basis of a first-principles-derived Hamitonian. ft is found that there exists a 
critical thickness of about three unit cells (~12 A) below which the ferroelectricity disappears under 
the condition that the in-plane polarizations are suppressed by sufficient clamping effect. Above the 
critical thickness, periodic 180° stripe domains with out-of-plane polarizations are formed in the 
systems in order to minimize the energy of the depolarizaing field. The stripe period increases with 
increasing film thickness. The microscopic mechanism responsible for these phenomena is discussed. 

PACS numbers: PACS:77.80.-e, 77.84.Dy, 77.80.Bh, 77.22.Ej 
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During the past decade, ferroelectric thin films have 
been investigated with a great deal of interest not only 
because their physical properties are different from those 
of bulk materials 1] but also because they are promising 
for microelectronic and micromechanical applications |2(. 
A fundamental problem attracting considerable atten- 
tions for ferroelectric thin film is the critical thickness for 
the ferroelectricity to occur. The observations that the 
ferroelectric ground states exist in 40-A-thick perovskite 
oxide film Qand in 10-A-thick crystalline copolymer Q 
seem to suggest the absence of the critical thickness. 
Some theoretical investigations supported that opinion 
by predicting ferroelectric ground states in various per- 
ovskite slabs |5|, Recently, Junquera and Ghosez in- 
vestigated a realistic ferroelectric-electrode interface and, 
contrary to the current thought, revealed a critical thick- 
ness of about six unit cells in BaTi0 3 thin filmQ. It 
was proposed that the depolarizing electrostatic field[j| 
is responsible for the disappearance of the ferroelectric- 
ity. However, in their study, only the monodomain case, 
where the atomic off-center displacements in every cell 
are the same, was considered. In fact, polydomains or 
other atomic off-center displacements may come into be- 
ing to remove the depolarizing electrostatic field. For 
example, Fu and Bellaiche found that a size reduction 
can lead to an unusual atomic off-center displacements 
vortex pattern in BaTiC>3 quantum dots Q. Therefore, 
it is worth investigating whether the film will form spe- 
cial atomic off-center displacements pattern to remove 
the depolarization field and whether the nonzero critical 
thickness still exists in that case. 

In this paper, by use of a first-principles-based ap- 
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proach, the ferroelectricity of disorder Pb(Zr . 5T10. 5)0-3 
thin films is investigated to clarify the role of the strain 
and the off-center displacement pattern in determining 
the critical thickness. It is shown that the ferroelectric- 
ity in thin films is closely related to the strain constraint 
imposed by the substrate. When the thin films are free 
from constraint, a nonzero polarization with an in-plane 
direction always exists in the system, which does not dis- 
appear even in the monolayer films. When there is a com- 
pressive strain, the in-plane polarization is suppressed 
while the out-of-plane polarization exhibits a strong de- 
pendence on the film thickness. Above a critical thick- 
ness of about three unit cells, the out-of-plane polariza- 
tion forms periodic 180° stripe domains, and the pe- 
riod increases with the film thickness. Below the crit- 
ical thickness, the stripe domain structure disappears, 
which suggests that the ferroelectricity can be suppressed 
even in the absence of the depolarizing electrostatic field. 
The orientation dependence of dipole-dipole interaction 
is revealed to be responsible for these phenomena, which 
is expected to generally work in perovskite ferroelectric 
films. 

We adopt the effective Hamiltonian of PZT alloys p ro- 
posed by Bellaiche, Garcia and Vanderbilt [101 fill fl3| to 
predict the properties of the disorder Pb(Zro. 5^0.5)03 
thin film surrounded by the vacuum. All the parameters 
of the Hamiltonian are derived from the first principle 
calculations and are listed in references [ToL ITl| . Cohen 
[l3| and Meyer et aL0] demonstrated that the effect 
of surface relaxation is significant in perovskites. How- 
ever, it should be noted that the effect of surface relax- 
ation has been considered implicitly in our simulations 
flEf , and the influence of the surface upon the ferroelec- 
tric order parameter is modest, as pointed out by Meyer 
et al. 14]. In our simulations, we do not include an 
external term of surface effect proposed by Fu and Bel- 
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laiche while simulating nanoscopic structures, since they 
demonstrated that the term has almost no effect on the 
polarization pattern^]- From the above effective Hamil- 
tonian, Monte Carlo simulations are conducted for large 
supercells typically containing between 5000 and 50000 
atoms mimicking the studied structures. The supercell 
average of the local soft modes (i is the cell index) 
is directly proportional to the macroscopic electrical po- 
larization. The influence of the substrate is imposed by 
confining the homogeneous in-plane strain. 

To efficiently calculate the long-range dipole-dipole in- 
teraction energy in thin films which lack the periodicity in 
the out-of-plane direction, we adopt the corrected three- 
dimensional Eward method, whose validity has been ver- 
ified analytically by Brodka and Grzybowski ^(| . In that 
scheme, a small empty space (about three times of the 
film thickness) introduced in the simulation box to sur- 
round the film can lead to very well converged results. It 
is noted that the Eward sum we adopted for the macro- 
scopic sphere is immersed in vacuum. Therefore the effect 
of depolarization field is considered in our model. 

In our simulation, the z axis ([001] direction) lies along 
the growth direction of the film, and the x and y axes are 
chosen to be along the pseudocubic [1 00] and [010] direc- 
tions. The film thickness d is measured by the number 
of unit cells along the z-axis. As the first step, we simu- 
late the ferroelectric properties of thin films without any 
constrain from the substrate, i.e., the free thin films are 
considered. It is found that for all the film thicknesses 
considered, there always exists a nonzero polarization 
along the x- or y- axis, and correspondingly, a tensile 
strain of about 2% along the same direction. (Data is 
not shown here.) It is consistent with the simulation re- 
sults in quantum dots 9] where the normal component of 
polarization near the surface is suppressed by the vac- 
uum. It also verifies the absence of any critical thickness 
for ferroelectricity in thin films under the condition of 
vanishing internal stress 0- 

In actual systems, the film strain will be imposed by 
the substrate (clamping effect) which has important ef- 
fects on the ferroelectricity [Pi ITsl Hij l2Cj . So we con- 
sider the strain effects in our simulations. We calculate 
the supercell average (u x , u y , u z ) of the local soft mode 
as functions of the in-plane strain. The result with d = 4 
is shown in Fig. I. It is clearly seen that u x and u y 
decrease with increasing (compressive) strain and disap- 
pear after the strain reaches a critical value, while u z 
keeps zero at any strain. The disappearance of the in- 
plane polarization is of the first order. It should be noted 
that when the in-plane strain is fixed, u x is equal to u y 
due to the symmetry requirement. Small difference be- 
tween u x and u v results from our treatment of avoid- 
ing effects of symmetry-equivalent rotations of the order 
parameter p^]. Detail simulations indicate that the crit- 
ical strain (~ 2%) is insensitive to the film thickness. 
From the above results, it can be concluded that the fer- 
roelectricity does not disappear in thin films (even for the 
monolayer case) if the strain imposed by the substrate is 
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FIG. 1: The dependence of the supercell average (u x ,u y , u z ) 
of the local soft mode on the strain for the film with the 
thickness of four unit cells (i.e., d — 4) 



not strong enough to suppress the in-plane polarization. 

Although the average value of the out-of-plane polar- 
ization (u z ) always keeps zero, the local mode u z exhibits 
completely different behavior below and above the criti- 
cal strain. u z is substantially small in any unit cell be- 
low the critical strain, while it quickly increases nearby 
the critical strain. Above the critical strain, periodic 
180° stripe domains are formed in the system. An ex- 
ample with stripe domains aligned along y-axis is shown 
in Fig. 2 for a thin film with d = 4 under 2% compres- 
sive strain. The domain walls are rather diffuse and the 
spatial variation of u z is close to a cosine function. The 
similar stripe domain stabilized by strain has also been 
found by Tinte et a2.|2l|. In the general cases, we find 
that the direction along which the stripe domains align is 
not certain. It may be the x or y axis and may also have 
an angle with the x or y axis. In each stripe domain, the 
z component (u z ) of local soft modes have the same sign. 
In other words, the stripe structure extends through the 
film in the z direction, which is important for ferroelectric 
field effect and has been observed in PZT[22j and PbTi03 
(PT) film • It is noted that the stripe domain only ap- 
pears above a certain film thickness. We calculate many 
kinds of supercells with d < 4. No domain structure or 
ferroelectricity is found in these supercells at 50 K even 
when the compressive strain is as large as 4%. Therefore, 
there is a critical thickness of about three unit cells (~1.2 
nm) in disorder Pb(Zro. 5^0.5)03 thin films below which 
the ferroelectricity disappears at large strains (> 2%). 

The dependence of ferroelectric-paraelectric phase 
transition temperature T c on the film thickness also sup- 
ports the existence of the critical thickness, as shown 
in Fig. 3. The temperature is rescaled as in Ref. [24j 
due to the fact that the effective-Hamiltonian approach 
overestimates T c . We find that T c , determined by the 
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FIG. 2: The y-z plane average (u x ,u y , and u z ) of the local 
soft mode as a function of x/a (a is the lattice constant) for 
120 x 5 x 4 supercell under 2% compressive strain. Only the 
data of 80 < x < 120 are shown. 



FIG. 3: Ferroelectric-paraelectric phase transition tempera- 
ture as the function of the film thickness d under 2% com- 
pressive strain. Dash line corresponds to T c of the bulk ma- 
terial. Inset: the polarization amplitude as a function of the 
temperature for d = 4. 



temperature dependence of the polarization amplitude 
(Fig. 3 inset), is far below T c of the bulk material (dash 
line in Fig. 3). This is consistent with the experimental 
observation of a significant decrease of T c in thinner films 
(below 20 nm)|23j. As the film thickness decreases, T c de- 
creases quickly. Extrapolation of T c in Fig. 3 shows that 
the ferroelectricity disappears at the thickness of about 
d = 3. It is interesting to note that Ghosez and Rabe 
also observed a critical thickness of 3 layers in PbTiC>3 
with completely different electrical boundary condition 
(perfectly screened condition) ,5] . Then, it seems that 
the electrical boundary conditions have little effect on 
the critical thickness, which is contrary to the common 
view. 

The stripe period for the film with d — 4 (Fig. 2) is 
determined to be 7.1 ± 0.5 cell units (~ 2.8 ± 0.2 nm) 
according to the fact that 17 periodic stripe domains are 
observed in an 120 x 5 x 4 supercell. We find that the pe- 
riod A of the strip domain has the same order of the film 
thickness d and increases with film thickness as shown in 
Fig. 4. 

Our predictions are helpful for understanding many 
experimental results. The stripe domains revealed here, 
with out-of-plane polarizations and 180° domain walls ex- 
tending through the film, are essentially identical to what 
have been observed in PT films by Streiffer et al. through 
the X-ray scattering measurements [23j . Their observed 
stripe period and the dependence on the film thickness 
are consistent with our calculated results. Based on our 
results, it is easy to understand why no detectable polar- 
ization was probed by the electrostatic force microscopy 
along the perpendicular directions of PZT filmj^, since 
it is demonstrated that the net polarization vanishes by 
forming the stripe domain with very short period. 

From a fundamental point of view, ferroelectricity in 
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FIG. 4: The effect of the film thickness on the stripe period 
under 2% compressive strain. The simulations for the period 
of stripe domain are conducted respectively for 120 x 5 x 4, 
160 x 5 x 6, and 180 x 5 x 8 supercells. The errors are also 
indicated in the figure. 



materials is determined by the competition between the 
long-range electrostatic and the short range covalent in- 
teractions. Our previous work revealed that the orien- 
tation dependence of dipole-dipole interaction is respon- 
sible for the appearance of unexpected phases in com- 
positional modulation PZT alloys near the morphotropic 
phase boundary j2f|. Here we indicate that the existence 
of the critical thickness and periodic stripe domains can 
be also understood by considering the orientation depen- 
dence of dipole-dipole interaction. The long-range polar- 
ization interactions are described in dipole-dipole inter- 
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action energy [Tllll2|: 



ing electrostatic field. However, whether the polarization 
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where Z* is the Born effective charge and the optical 
dielectric constant of the material. To simplify the the- 
oretical analyzes, we assume that all are along the z 
axis. For the monolayer film, the dipole-dipole interac- 
tions are constrained to be intralayered. Then the second 
term of Eq. (1) is zero and the dipole-dipole interaction 
coefficient is always positive for the ferroelectric phase. 
Therefore, in the monolayer film, the dipole-dipole inter- 
action does not support the formation of ferroelectricity 
along the z axis. The ferroelectric phase will not appear 
in this case unless the in-plane polarization is considered, 
as what we have demonstrated above. 

In contrary to the intralayer dipole interaction, the in- 
terlayer dipole interaction is in favor of the ferroelectric 
order, and its effect becomes more and more important 
with increasing film thickness. So after the film thick- 
ness reaches a certain critical value, the ferroelectricity 
appears. If the depolarization field has not been elim- 
inated by such effects as the surface charge screening, 
domain structures will occur to decrease the intralayer 
dipole interaction energy. 

The variation of the stripe period with the film thick- 
ness comes from the delicate balance between the long- 
range dipole interaction and the short-range interaction. 
The dipole-dipole interaction energy increases with in- 
creasing modulation period, while the short-range inter- 
action energy decreases with increasing period. Com- 
pared with the short-range interaction, the dipole inter- 
action is long-range and thus more closely related to the 
film thickness. This feature determines that the stripe 
period varies with the film thickness. To further clarify 
this point, we present a simple numerical analysis in the 
following. The local modes are assumed to vary along 
the x axis as 



u z0 cos( 
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where A is the period of the strip domain and u z o is the 
amplitude of local mode u z . In Fig. 5, the average dipole- 
dipole interaction energy and short-long interaction en- 
ergy of one unit cell are shown as functions of A for thin 
films with different thickness. The larger the period is, 
the more quickly the dipole-dipole interaction energy de- 
creases with increasing film thickness. On the contrary, 
the energy difference between different modulation peri- 
ods is almost independent of the film thickness for the 
short-range interaction, i.e., increasing the film thickness 
only shifts the curve upward. Therefore, the minimum of 
the total energy (see the inset in Fig. 5) moves to higher 
A values while increasing the film thickness. 

The formation of domain structure helps to decrease 
the dipole interaction energy by screening the depolariz- 
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FIG. 5: The dipole-dipole interaction energy (dash line) and 
short-range interaction energy (solid line) as functions of the 
stripe period A for the thin film with different thickness d. 
Inset: The sum of the dipole-dipole interaction energy and 
short-range interaction energy as functions of the stripe pe- 
riod for d — 4. 



will appear is also affected by the short-range interaction. 
If the minimum of the total energy of domain structure 
is higher than that of the paraelectric phase, the fer- 
roelectricity will be suppressed, which is just the cases 
occurring below the critical film thickness (note that the 
total energy shown in Fig. 5 will increase with decreasing 
thickness). It suggests that the critical thickness exists 
even if the 180° domains are formed to screen the depo- 
larizing field. 

In summary, we have demonstrated the existence of 
periodic 180° stripe domains with the out-of-plane po- 
larization and a critical thickness for ferroelectricity in 
disorder PZT thin films, which may widely exist in per- 
ovskitc thin films. It seems that the electrical boundary 
conditions have little effect on the critical thickness since 
the critical thickness is almost the same for completely 
different electrical boundary condition. It is pointed out 
that the critical thickness appears only when the strain 
imposed by the substrate is strong enough to suppress 
the in-plane polarization. Otherwise, ferroelectricity can 
be maintained even in monolayer film. 
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